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ABSTRACT 

Results are presented for [CII] 158 /im line fluxes observed with the Herschel 
PACS instrument in 112 sources with both starburst and AGN classifications, 
of which 102 sources have confident detections. Results are compared with mid- 
infrared spectra from the Spitzer Infrared Spectrometer and with L ir from IRAS 
fluxes; AGN/starburst classifications are determined from equivalent width of the 
6.2 ixm PAH feature. It is found that the [CII] line flux correlates closely with the 
flux of the 11.3 ^m PAH feature independent of AGN/starburst classification, 
log [f([CII] 158 jan)/f(11.3 ^m PAH)] = -0.22 ± 0.25. It is concluded that 
[CII] line flux measures the photodissociation region associated with starbursts 
in the same fashion as the PAH feature. A calibration of star formation rate for 
the starburst component in any source having [CII] is derived comparing [CII] 
luminosity L([CII]) to L ir with the result that log SFR = log LQCII)]) - 7.08 ± 
0.3, for SFR in M yr" 1 and LQCII]) in L . The decreasing ratio of LQCII]) to 
L ir in more luminous sources (the "[CII] deficit") is shown to be a consequence 
of the dominant contribution to Li T arising from a luminous AGN component 
because the sources with largest L ir and smallest LQCII] )/L ir are AGN. 
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Introduction 



Understanding the initial formation of galaxies depends on discovering sources obscured 
by dust and tracing these sources to their earliest epoch in the universe. The extreme 
luminosity of dusty, local sources was originally revealed by the Ultraluminous Infrared 
Galaxies (ULIRGs, e.g. Soifer, Neugebauer and Houck 1987, Sanders and Mirabel 1996), 
whose luminosity arises from infrared emission by dust, and this dust often obscures the 
primary optical sources of luminosity. That such galaxies are important in the early universe 
was demonstrated by source modeling which indicate d that the infrared dust emission from 



galaxies dominates the cosm ic background luminosity (jChary and Elbazll200ll ; lLagache et al. 
20041 ; iLe Floc'h et alJliooih . 



Surveys in the submillimeter were the first to disc over individual, optically obscured, 
dusty sources at redshifts z > 2 ([Chapman et al.ll2005l ). A variety of observing programs 



quently found luminous ULIRGS to redshifts z ~ 3 (e.g. 



Houck et al. 


2005; 


Yan et al. 


2007 



Sajina et al.ll2007l ; IWeedman and Houckl|2009bl ). This S pit zer- discovered population of high 



redshift ULIRGs has large infrared to optical flux ratios [/^(24 /im) > 1 mJy and R > 24] at- 
tributed to heavy extinction by dust and has been labeled "dust obscured galaxies" (DOGS; 
Dey et al. 2008). Some DOGs are powered primarily by starbursts and some by active galac- 
tic nuclei (AGN), and the DOGS are similar to the population of submill imeter galaxies in 



overall spectral energy distribut i ons (SEDs), redshi f ts, and luminositie s (IPope et al 
Menendez-Delmestre et al.ll2009i ; ICoppin et al.ll2010i ; iKovacs et al.ll2010l ). 



2008 



To discover and understand dusty galaxies at even higher redshifts than the DOGs 
known so far, the atomic line emission of [CII] 158 /im is t he single most important spectro- 



scopi c feature because it is the strongest far- infrared line (jStacey et al.lll99ll ; iLuhman et al. 



20031 ; iBrauher et al.l 120081 ) . As a consequence, this line will provide the best opportu- 
nity for redshift determinations and source diagnostics using submillimeter and millime- 
ter spectroscopic observations. Already, [C II] has been detected at redshift exceeding 7 



(IMaiolino et al. 



2005; 



z < 2.5 Irlailev-Dunsheath et aDboiol : Istacev et al.l fioiol : llvison etaDboid ). 



Venemans et al. 12011) and shown to be strong in starb ursts with 1 < 



Our primary motives for the present paper are to present [CII] results for a large sample 
of dusty sources and to compare with mid-infrared classification indicators for starbursts 
and AGN. This comparison leads to a calibration between star formation rate (SFR) and 



1 Based on observations with the Herschel Space Observatory, which is an ESA space observatory with 
science instruments provided by European-led Principal Investigator consortia and with important partici- 
pation from NASA. 
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[CII] luminosity. We emphasize the diagnostics used for DOGS at z ~ 2, because the large 
populations of submillimeter and mid-infrared DOGS now known at this epoch provide a 
crucial reference for scaling to higher redshifts. The epoch 2 < z < 3 is also important 



because this is the observed epoch at which starburst and AGN activity seems to peak (e.g. 


Madau et al. 


1998; 


Reddv and Steidel 


2009; 


Fan et al. 


2004; 


Croom et al. 


2004; 


Brown et al. 



2006|). 



The [CII] line should be primarily a diagnostic of star for mation, being associated with 



the photodissocia t ion region (PDR) su r rounding starbursts ( Tielens and Hollenbachl 11985 



Helou et al.l l200ll ; iMalhotra et al.l l2001t iMeijerink et al.l 120071 ). and the line appears to be 
weaker in the most luminous sources ("the [CII] deficit", Luhman et al. 2003). It is crucial 
to understand the origin of this line and the extent to which its luminosity is a measure of 
SFR. Does [CII] scale with other star formation indicators? Is the [CII] deficit a consequence 
of AGN dominance rather than star formation in luminous sources? Determining such 
answers is th e objective of new o bservations we have undertaken with the Herschel Space 
Observatory Jpilbratt et al.lboioh . 



2. Observations 



2.1. Selection of Sources 



Using the Herschel PACS instrument (IPoglitsch et al.ll2010l ). it is now possible to mea- 
sure efficiently the [CII] luminosity in luminous, dusty galaxies within the local universe. 
Many observations are underway (e.g. the SHINING key program, E. Sturm, P.I.). Our 
Herschel PACS observing program (OTldweedmanl) includes 112 sources chosen to con- 
nect the [CII] results to various mid-infrared diagnostics of starburst and AGN activity that 
can be derived from spectroscopic observations with the Spitzer IRS. 

Our source list was assembled using these criteria: 1. All targets have complete low 
resolution and high resolution mid-infrared spectra from 5 /xm to 35 /im with the Spitzer IRS; 
low resolution spectra are available in the Cornell Atlas of Spitzer IRS Spectra (CASSIS; 
Lebouteiller et al. 20110). 2. All targets have complete fluxes from the Infrared Astronomical 
S atellite (IRAS) so the total infrared luminosities Lj r can be determined with the relation 
of ISanders and Mirabel! (119961 ). 3. Finally, all targets are spatially unresolved (according to 
estimates described below) to give confidence that IRS and Herschel spectroscopy measure 



2 http://cassis. astro. cornell.edu/atlas; CASSIS is a product of the Infrared Science Center at Cornell 
University. 
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the same source. These selection criteria do not include any source classification criteria; the 
criteria derive only from selection based on available archival observations. 



The 112 sources in our observing program are taken from the 301 sources in lSargsyan et al 



( 120111 ) by proceeding as follows. That list was produced starting with all IRS archival obser- 
vations then available having both IRS low resolution spectra and complete flux measures 
with IRAS, giving a sample of 501 sources. To exclude extended sources for which IRS and 
IRAS flux comparisons would not be the same, the total f v (IRAS 25 fim) was compared to 
the fj,(IRS 25 fim) measured with the 10" slit of the IRS t o def ine a list of 301 sources 



estimated to be unresolved (see Figure 2 in ISargsyan et al.l (120111 ) and the accompanying 



discussion). The resulting sample has 0.004 < z < 0.34 and 42.5 < log Ljr < 46.8 (erg s x ) 
and covers the full range of starburst galaxy and AGN classifications. 

Of these 301 sources, 182 also have IRS high resolution spectra in addition to the low 
resolution spectra. At the time of our PACS proposal, 41 of these 182 had PACS observations 
listed either in the Herschel Reserved Observations Search Tool (30 sources) or archival 
[CII] data from ISO (11 sources; Brauher et al. 2008). Of the 141 remaining sources, we 
selected for new PACS [CII] observations the 123 brightest as measured by IRS flux of the 
11.3 fim polycyclic aromatic hydrocarbon (PAH) feature. Subsequently, 11 of these 123 
were yielded to program OTldfarrahl, giving our final sample of 112 sources. These are 
the results we report in the present paper. We also in clude for comparison the prototype 



ULI RG Markarian 231 usi ng the PACS [CII] flux from Fischer et all feOlOf ) and IRS data 



from lSargsyan et al.l (120111 ) because this important source satisfies all of our sample selection 



criteria. 

Having IRS low resolution and high resolution spectra for all sources means that many 
comparisons can be made with various mid-infrared diagnostic features, including several 
atomic and molecular emission lines, silicate absorption and emission, and the dust con- 
tinuum. In future papers, we will present comparisons with emission line fluxes, velocities 
and profiles. For the present analysis, our goal is only to compare the [CII] results to PAH 
molecular emission features. The PAH features are the most important diagnostics for clas- 
sifying dusty, luminous starbursts and AGN at z > 2, as described below in section 2.2 and 
illustrated in Figure 1. They are also the most important mid- infrared measure of SFR. By 
relating [CII] to PAH, the populations of DOGS already known at z ~ 2 can be compared 
with sources at even higher redshifts when such sources are observed using [CII]. 
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2.2. Classification of Sources 



For understanding the ultimate luminosity source, it is necessary to distinguish star- 
bursts and AGN. The most important mid-infrared spectroscopic criterion for classification 
is the broa d complex of PAH e mission (Figure 1) which arises within the PDR surrounding 
starbursts ( IPeeters et al.ll2004l ). That PAH em ission is weak i n AG N compared to starbursts 
was initially demonstrated observationally by iGenzel et al.l (119981 ) using spectroscopy with 
the Infrared Space Observatory (ISO). The strength of a PAH feature compared to the un- 
derlying continuum (the equivalent width, EW) decreases as the AGN component increases 
because the continuum beneath the PAH feature increases in proportion to the hot dust 
heated by the AGN. 

Spectra for hundreds of sources with the IRS quantified the depende nce of PAH strength 



on starburst/AGN classifica t ion determined from optical spectra (e.g. iBrandl et al.l 12006 
Weedman and Houckl l2009al ; ISargsyan et al.l 1201 ll ). These comparisons with optical classi- 
fications led to an infrared classification derived from EW(6.2 /im), measured as described 
below. The continuous gradation of EW(6.2 /iin) provides a single parameter, quantitative 
measure of the starburst/AGN mix. We use the EW(6.2 /im) for classification rather than 
the stronger 11.3 /iin PAH feature because the 11.3 /im feature is near the long wavelength 
limit of the IRS for z > 2, and we also want to apply the classification at high redshifts. 

The PAH spectroscopic emission features are complex (Figure 1), including an under- 
lying plateau between 5 /im and 10 //m, a maximum at rest-frame ~ 7.7 /im , and specific 
features at 6.2 /im and 11.3 /im. Features are broad and blended, and the "continuum" 
beneath the features includes wings of other PAH emission features. Sophisticated measures 
of PAH strengths require various assumptions about the underlying dust continuum and rel- 
ative feature stren gths to deconvolve the full, broad PAH features into different components 
(jSmith et al.ll2007l ). This deconvolution is uncertain when using observations of faint sources 
with poor signal to noise ratios (S/N) and limited coverage of rest frame wavelengths, as 
arises for the highest redshift sources with IRS spectra. 

To avoid such uncertainties, we use simple parameters for PAH strength that allow 
consistent, model-independent observational measures that are possible even with weak fea- 
tures or poor S/N. These measures are made with the SMART software for IRS spectra 
(IHigdon et al.l 120041 ) and are: 1. the EW(6.2 /im) determined using a gaussian fit to the 
PAH feature and a linear fit to the continuum beneath the feature within the range 5.5 
/im to 6.9 /im; 2. the total flux of the 11.3 /im feature, f(11.3 /im), determined with a 
gaussian fit to the PAH feature and a linear fi t to the continuum be neath the feature from 
10.5 /im to 12 /xm. Measurements are given in ISargsyan et al.l ( 1201 if ). 
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4 6 8 10 12 14 

rest wavelength (/^m) 

Fig. 1. — Transition in mid-infrared spectra from pure starburst to pure absorbed AGN. 
Top spectrum is median observed rest frame IRS spectrum of 51 starburst PAH emission 
sources (100% Starburst) from samples in Sargsyan et al. (2011); bottom spectrum is median 
observed rest frame IRS spectrum of 65 AGN with silicate absorption (100% AGN; from same 
reference); intermediate spectra show mixes of top and bottom spectra, changing mix by 25% 
in each spectrum. Spectra are normalized to peak f u (7.8 //m) and displaced by 0.5 units of 
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fraction of luminosity from AGN 

Fig. 2. — Equivalent width in fj,m of 6.2 fj,m PAH emission feature shown in Figure 1 
compared to fraction of total luminosity at 7.8 fim which arises from AGN component. Curve 
derives from various mixtures of the "100% Starburst" and "100% AGN" spectra shown in 
Figure 1. Criteria for EW(6.2 ^m) used to classify Herschel sources in the remaining figures 
are illustrated. 
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Although the range of EW(6.2 /xm) is continuous among spectra, comparisons with 
optical classes show consistent divisions. The initial study of IRS spectra for 22 opti- 
cally classifi^d_j^rbmsts showed that 21 of 22 sources have rest frame EW(6.2 /im) > 0.4 



/im (IBrandl et al.ll2006l ). Subsequent studies confirmed this result and also determined an 



EW limit for optically classified AGN. Of the 19 sources with EW(6.2 /xm) > 0.4 /xm havin g 



optical classifications in the f„(24 /xm) > 10 mJy sample of IWeedman and Houckl (j2009al ) 



18 are classified as starbursts ; of the 17 sources wit h EW(6.2 /xm) < 0.1 /im, 16 are classified 



as AGN. The distribution in lSargsyan et al.l (120111 ) is illustrated in their Figure 4, where 14 



of 17 sources having both EW(6.2 /im) > 0.4 /im and optical classifications are starbursts, 
and 42 of 57 having EW(6.2 /xm ) < 0.1 /xm are AGN. Using these results, we have adopted 
the criteria that AGN have EW(6.2 /xm) < 0.1 /xm, composite sources have intermediate 
0.1 /im < EW(6.2 /xm) < 0.4 /xm, and starbursts have EW(6.2 /xm) > 0.4 /xm. (Observed 
frame equivalent widths are greater by a factor of (1+z) compared to rest frame equivalent 
widths.) 

This empirical observational classification is quantitatively consistent with starburst/AGN 
luminosity fractions when synthetic spectral combinations of mid-infrared spectral templates 
are assumed. Using the mixes of observed starburst and AGN spectra shown in Figure 1, 
for example, the dependence of EW(6.2 /xm) on the fractional AGN contribution to the 
mid-infrared luminosity is shown in Figure 2. With this mix of spectra, > 90% of the mid- 
infrared luminosity arises from a starburst if EW(6.2 /xm) > 0.4 /xm, confirming that this 
value of EW(6.2 /xm) defines a starburst. The majority of the mid-infrared luminosity, > 
55%, arises from an AGN if EW(6.2 /xm) < 0.1 /im . In fact, most AGN in the present paper 
actually have EW(6.2 /im) < 0.01 /xm for which > 90% of the luminosity arises from an 
AGN in Figure 2. Composite sources have contributions from both starburst and AGN. The 
boundaries illustrated in Figure 2, therefore, are those used to define classification symbols 
in those figures below which do not display a quantitative EW(6.2 /xm). 



2.3. Observations and [CII] Measurements 

All [CII] observations were made using PACS line spectroscopy in point source chop nod 
mode with medium throw. A single repetition cycle was used for all but 20 sources, giving 
total observing time for the program of 20.2 hours for 112 sources. Table 1 includes results 
for the 112 sources. Data reduction was done with version 8 of the Herschel Interactive 



Proc essing Environment (HIPE), together with the "PACSman" software (ILebouteiller et al. 



20121 ) used for fitting line profiles and continuum within each of the 25 equivalent spatial 
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pixels, or "spaxels", produced by the PACS image sliceiQ Illustrations of these fits are in 
Figure 3, with examples of both high and low S/N sources, together with an example of a 
source with an unusually broad line profile. 

We include a line detection as real if the S/N for the line profile fit in the brightest 
spaxel exceeds 3. If S/N < 3 for the line flux in the brightest spaxel, we adopt an upper 
limit for the total line flux. Of our 112 sources, 102 are detected according to this criterion. 
All results are in Table 1. 

One method for measuring total source fluxes with PACS is to take only fluxes in 
the brightest, central spaxel and correct to total flux assuming an unresolved source that 
is perfectly centered in the central spaxel. In this case, the flux from the central spaxel is 
corrected by a factor of 1.95 at wavelength 160 |im to deter™ the total flvefl The greatest 
uncertainties in this technique are the requirement of perfect centering and the possibility 
that sources are extended. 

To minimize centering uncertainties and to maximize the inclusion of total flux, the 
procedure we adopt instead for measuring detections is to measure the total line flux in the 
brightest spaxel plus the 8 surrounding spaxels, f(3x3). (For detected sources, we find that 
the brightest spaxel is always the central spaxel, except for 4 cases where it is displaced by 
one spaxel, noted in Table 1). The f(3x3) covers a spatial area of 28" by 28". Calibration 
communicated to us by the PACS calibration team gives the result that the f(3x3) should 
be corrected by a factor 1.16 to give the total flux at 160 /im for an unresolved source 
that is precisely centered. This correction is slightly wavelength dependent, increasing to 
a maximum correction of 1.21 for our highest redshift source, for which the [CII] line is 
observed at 212 fim. 

To list the total fluxes of detected sources, therefore, we correct f(3x3) to include flux 
outside these spaxels by taking values between 1.16f(3x3) and 1.21f(3x3) as the total source 
line flux, depending on the observed frame wavelength according to the formula given in the 
footnote to Table 1. These corrected f(3x3) are the values listed in Table 1 for the total [CII] 
line flux. 

These two alternative measures of total line fluxes from the brightest spaxel and from 
f(3x3) are compared in Figure 4. This comparison illustrates empirically the large error that 
can arise from imperfect centering if using only the brightest spaxel. The plot compares the 



3 http: / /herschel.esac.esa.int /Docs /PACS /pdf/pacs-om.pdf 

4 http: / /herschel.esac.esa.int / twiki / pub/Public/PacsCalibrationWeb /PacsSpectroscopyPerformanceAndCalibration- 
31May2011.pdf 
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Fig. 3. — Examples of line fits in the brightest spaxel for three sources in Table 1. Vertical 
axes are flux density in Jy; horizontal axes are wavelength differences in //m from [CII] rest 
wavelengths derived using redshifts in Table 1. Upper left panel (source 28 in Table 1) is 
a high S/N detection for which line flux derives from gaussian fit shown; upper right panel 
(source 92 in Table 1) is an upper limit; lower panel (source 83) is a complex, broad profile 
for which line flux is total integrated flux within the wavelength bounds of the two profiles 
shown. [CII] line fluxes in Table 1 are the sum of the brightest spaxel plus the 8 surrounding 
spaxels, increased by a correction factor of 1.16 to 1.21 (depending on redshift) to include 
the flux from an unresolved source that would fall outside these 9 spaxels. 
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Fig. 4. — Vertical axis compares total corrected [CII] source flux derived from observed 
flux in 3x3 spaxels compared to corrected flux derived only from the brightest spaxel; both 
observed measures are corrected for fractional flux that would fall outside of those spaxels for 
a perfectly centered, unresolved source, using corrections of 1.95f (brightest) and 1.16f(3x3) 
at 158 jum (small changes in corrections because of different redshifts are ignored in this 
ratio comparison). Horizontal axis is corrected total [CII] flux 1.16f(3x3) in W m~ 2 . The 
102 detected sources in Table 1 are shown as crosses for AGN, open circles for composite 
AGN plus starburst, and asterisks for starbursts (all using the EW(6.2 /im) classification 
discussed in section 2.1). Vertical line defines flux upper limit for the 10 undetected sources 
in Table 1. Upper horizontal line is maximum value of ratio 1.16f(3x3)/1.95f(brightest) that 
can arise for an unresolved source because of imperfect centering; sources above this line 
must be extended. Lower horizontal line is ratio of unity expected for a perfectly centered 
unresolved source. 
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single brightest spaxel corrected by a factor of 1.95 with the sum of the 9 spaxels corrected 
by 1.16. The limiting ratio for these corrected fluxes should be unity (log = on vertical 
axis) for a perfectly centered point source, and several sources satisfy this value. Small 
centering uncertainties have little effect within the large area of the f(3x3) flux, so the 
scatter in 1.16f(3x3)/1.95f(brightest) primarily demonstrates the uncertainty arising when 
using only the brightest spaxel to determine total flux. The vertical dispersion of the points 
and the dominance of values above unity can be explained by slight differences in source 
centering within the brightest spaxel. (To simplify this illustrative calculation, we apply no 
wavelength dependence for different redshifts to the correction factors for Figure 4 because 
such differences are small compared to the centering uncertainties.) 

For example, in the extreme case of a source being offset 5" in a direction such that 
the centering is on a spaxel corner, the "brightest spaxel" would be shared evenly among 
4 spaxels. In this extreme example, the brightest spaxel would have 1/4 of the flux within 
the 4 spaxels, or within encircled energy of radius 10" (taken as 72% of the total flux for 
an unresolved source), so the brightest spaxel has 18% of the total source flux. For this 
extreme case of sharing the brightest spaxel among 4 spaxels, therefore, the correction to 
total flux derived only from the brightest spaxel should actually be a factor of 5.6 and not 
a factor of 1.95, so the corrected total flux from the brightest spaxel alone if corrected by 
1.95f(brightest) would be erroneously faint by a factor of 2.9. This extreme case is shown 
as the upper line in Figure 4. Intermediate centering of an unresolved source within the 
brightest spaxel would lead to values anywhere between the upper and lower lines in Figure 
4, so centering uncertainties alone of < ± 5" could explain the dispersion of all points in 
Figure 4, except for the one source falling above the upper line. 

Extended sources would also have ratios above unity in Figure 4, and it is not possible 
to determine from this ratio alone whether a source is unresolved but not perfectly centered, 
or whether the source is slightly extended. A diffuse source so extended that it evenly fills 
the central 3x3 spaxels (so each of these 9 spaxels has the same flux as the brightest spaxel) 
would have log[1.16f(3x3)/1.95f(brightest)] = 0.73, but no sources show this extreme ratio. 
This is proof that no sources are extended as large as 30". 

Further evidence that most sources are unresolved is the similar distribution in ratio 
1.16f(3x3)/1.95f(brightest) between starbursts and AGN. If some sources are extended, they 
should be extended starbursts instead of unresolved AGN, but starbursts show no more 
extension to large ratios than do the AGN. The consequence to our analysis of having some 
marginally resolved sources would be that the total [CII] flux is erroneously large compared 
to PAH, because the PAH measure derives from IRS measures calibrated for unresolved 
sources. When comparing the [CII]/PAH ratio below in section 3.2, however, the only 3 
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sources with [CII]/PAH ratios significantly above the dispersion are three AGN, and none 
of these three show any evidence of spatial extent in the PACS images. 

Figure 4 also illustrates how upper limits for non detections are determined. The detec- 
tion limit for our sample is at log f( [CII] 158 /im ) ~ -16.7 in units of W m -2 . This value is 
taken as the upper limit of [CII] line flux for the 10 sources which are not detected, according 
to our criterion of S/N < 3 in the brightest spaxel. 



3. Analysis and Results 



In the present analysis, the mid-infrared spectroscopic features used for comparisons to 
[CII] line luminos ities are the PAH emis sion features at 6.2 /im and 11.3 /im . These values 
are tabulated by ISargsyan et al.l (120111 ) measured as described above in section 2.1. The 
11.3 /im PAH is used as the PAH flux measurement for comparison with [CII] because this 
feature is detected in all but two of our sources, whereas the 6.2 /im feature is only a limit 
in 25 sources. The EW(6.2 /im) is used only for source classification, and all sources with 
limits for EW(6.2 /im) have classification as AGN because of the weak 6.2 /im PAH. 



3.1. Comparison of PAH and [CII] 



Photoionization and photoelectric heat ing models for [CII] (jStacey et al.lll99ll ; iHollenbach and Tielens 



1999; Malhotra et all 



tions (ILuhman et al 



200 ll ; Abel et al. 2009^ in comparison with previousl y available observa 



20031 : iGracia-Carpio et al.ll201ll ; IStacev et aUl2010h indicate that [CII] 



emission arises primarily in PDRs, with ionization produced by the hot stars of the adjacent 
HII region. Although [CII] can arise in any region with singly ionized carbon and sufficiently 
energetic electrons for collisional excitation, the level of ionization seems the most important 
parameter for line strength; models show a weakening of [CII] if the ionization parameter 
increases, resulting in a diminishing of the PDR compared to the HII region. This could 
be the result of harder ionization either from AGN or from unusually hot stars in compact, 
young starbursts. 

Comparisons of 6.2 /im PAH and [CII] using ISO results showed overall consistency 
between the two measures of PDRs but discussed vari ous reasons why detailed agreement 
is not expected (IHelou et al.ll200ll ; ILuhman et al.l 120031 ). Models as well as Herschel obser- 
vations of spatially re solved [CII] emission regions and Galactic PDRs show that detailed 



relations are comp lex (IMookerjea et al.ll2011l ; iLebouteiller et al.ll2012t iVelusamy et al.ll2012 
Beirao et al.ll2012l ) but generally confirm that [CII] emission arises primarily in the PDR and 
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scales with the PAH flCroxall et all GO 12 ). 



For our sample, the objective is an empirical observational comparison to determine 
how the observed [CII] and PAH fluxes scale together when integrated over many star form- 
ing regions throughout many sources. Such scaling would be an indication that the [CII] 
luminosity can be used as a quan titative measure of PDRs and SFR in the same fashion as 
the PAH luminosity can be used (IPeeters et al.ll2004l ) and thereby allow reliable use of [CII] 
as a star formation indicator. The large sample of sources enabled by the new Herschel 
observations can also determine the cosmic dispersion in this comparison and the extent to 
which the differen ces in [CII] /PAH ratios seen on small scales within starburst galaxies (e.g. 
Beirao et al.ll2012[ ) average out when integrated over entire starburst systems. 



Figure 5 compares the [CII] total line flux, f( [CII] 158 /mi), with the 11.3 fim PAH 
feature, f(11.3 /zm PAH), as a function of the source classification from EW(6.2 /an). Figure 
6 shows the [CII] /PAH ratio compared to total infrared luminosity L,>. (For comparison, the 
groun d based measurement of [CII] for the source at z = 1.3 from lHailey-Dunsheath et al. 



( 120101 ) is also shown in Figure 5 because this source has an IRS spectrum allowing PAH 
measurement; this source has higher redshift than any of our PACS sources.) Also shown 
for comparison in thes e and subsequent figu res is the prototype local ULIRG Markar i an 23 1 
using [CII] fluxes from lFischer et all ( 120101 ) and PAH measures from ISargsyan et all ( 1201 ll ). 
because this source satisfies our selection criteria. 

Figure 5 shows the important conclusion that the [CII] to PAH ratio is independent of 
starburst/AGN classification. The median ratios for all classes are the same. Figure 6 shows 
that the ratio does not depend on source luminosity. The overall median and dispersion 
for all sources is log [f([CII] 158 /W)/f(11.3 /mi PAH)] = -0.22 ± 0.25. This scatter is a 
measure of the cosmic dispersion in the comparison of [CII] and PAH when integrated fluxes 
are observed that include entire, unresolved collections of starbursts. This dispersion is 
independent of AGN/starburst classification which implies that the dispersion is a measure 
of variances within the starbursts. 

Because the PAH luminosity is determined by the starburst, this result means that the 
[CII] line flux also depends only on the starburst component of the source, regardless of the 
fractional starburst/AGN mix. This result indicates that [CII] luminosity can be used as a 
measure of starburst luminosity with as much reliability as PAH, because [CII] is a measure 
of the same PDRs as is the mid-infrared PAH. The [CII] luminosity measures SFR in any 
source for which [CII] is observed, even if the source luminosity is dominated by an AGN. 

A question about use of the 11.3 /im PAH feature for flux measures is whether this 
feature might be affected by silicate extinction in sources showing strong silicate absorption 
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Fig. 5. — Ratio of [CII] to PAH 11.3 /im line fluxes, compared to source classification from 
EW(PAH 6.2 ^m) measured in jj,m. Crosses are AGN from the EW classification (discussed 
in section 2.1); thick crosses are those sources with silicate absorption in IRS spectra noted 
in Sargsyan et al. (2011). Open circles are composite AGN plus starburst, and asterisks are 
starbursts. Sources with diamonds (all AGN) are upper limits to [CII] line fluxes in Table 1. 
Triangle is the source at z = 1.3 from Hailey-Dunsheath et al. (2010). Square is Markarian 
231 using [CII] flux from Fischer et al. (2010). Horizontal bars are medians within each 
category; medians include limits because all limits fall below the median. Vertical error bar 
shows the observational uncertainty for individual points assuming flux uncertainties of ± 
12% for f([CII]) and ± 10% for f(PAH 11.3 /im). 



- 16 - 



1 .0 P 



< 

CL 



=1 
CO 



cn 

o 



0.5 



0.0 



-0.5 



9 



o 

+ + + + ^ 
o o O 

° + + 9^ O * 

o 



o 



O 



' + 



+ o o + 
"o 



1 



11 12 
log L(IR; 



1 3 



1 4 



Fig. 6. — Ratio of [CII] and PAH 11.3 /im line fluxes, compared to total infrared luminosity 
L ir in L . Crosses are AGN from the EW classification, open circles are composite AGN 
plus starburst, and asterisks are starbursts. Sources with diamonds (all AGN) are upper 
limits to [CII] line fluxes in Table I. Square is Markarian 231. 



-17- 



as in Figure 1, because the PAH feature falls within the broad silicate absorption centered 
at 9.7 /mi. The silicate absorption is always associ ated with AGN clas s ification and is 



attrib uted to thick dust clouds close to the AGN (e.g. Ilmanishi et al .1120071 ; iLevenson et al. 



20071 ). An important question about the geometry of these absorbing clouds is whether they 
also surround the starburst regions. The [CII] results allow a test of this. If the starbursts are 
affected by the same silicate absorption that affects the AGN, then the f( [CII] 158 /xm)/f(11.3 
jum PAH) ratio should be systematically larger in sources with silicate absorption because 
the [CII] is not affected by the silicate absorption feature. 

These absorbed sources are noted in Table 1 and shown in Figure 5, where it is seen that 
there is no systematic difference in the ratio for the absorbed AGN compared to the emission 
AGN. Formally, the median for the absorbed AGN is log [f([CII] 158 /W)/f(11.3 fjem PAH)] 
= -0.23 and is -0.18 for the emission AGN, which is in the opposite sense that would arise if 
the absorbed AGN have suppressed 11.3 /im PAH. (The AGN median shown in the Figure 
is for all AGN.) Based on this observation, we conclude that the 11.3 fim PAH feature arises 
outside of the region affected by silicate absorption, so this feature has the same reliability 
for measuring starburst PDRs regardless of the presence of silicate absorption. 

Among all of the sources in Figures 5 and 6, only three appear anomalous. These are all 
AGN and have the largest ratios f([CII])/f(PAH) - sources 3C 120, Mrk 590, and NGC 3393 
in Table 1 - with f([CII])/f(PAH) about a factor of ten greater than the median. All have 
strong and reliable [CII] detections with no unusual profiles or evidence of spatial mismatches 
or source extension. The PAH 11.3 /im feature is weak in all three, but independent CASSIS 
spectra yield fluxes consistent to within ~ 30%, so there is no indication that the anomalous 
weakness of PAH is a measurement error. All are silicate emission sources, so the anomalies 
cannot be explained by silicate absorption associated with the AGN. These are the only 
sources among the 112 sources in all categories with excess [CII] luminosity compared to 
PAH. Why are the ratios so unusual in these three sources? 

There are various possible explanations. One possibility is that an unusual combination 
of radiation hardness and ionization parameter, caused by geometry, density distributions 
or ionizing spectrum, caus es PAH emission in t he PDR to be suppressed while maintaining 



[CII] within the HII region ( ILuhman et al.ll2003f ). An alternative possibility could be that the 



PAH are weak because of star formation taking place in dense clouds which are so obscured 
that the mid-infrared PAH from the starburst PDR suffers significant extinction compared 
to the far infrared [CII]. Because the AGN is not absorbed, such clouds would require a small 
filling factor or geometric placement outside of the observer's line of sight to the AGN. In 
this circumstance, we would expect mid-infrared emission lines associated with any obscured 
starburst, such as [Ne II] 12.8 /mi, also to be unusually weak compared to [CII], and would 
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also expect to observe excess far infrared continuum from the obscured starbursts. We defer 
further analysis of this question until making more comparisons among Spitzer mid-infrared 
emission lines, [CII] and SEDs for our full sample of sources. 



3.2. Star Formation Rate from [CII] 



For eventual application to high redshift, dusty sources in which [CII] measurements 
with ALMA or other submillimeter /millimeter spectroscopy are the primary diagnostic, cal- 
ibration of [CII] luminosity to star formation rate is our most important objective. The 
conclusion reached above, that [CII] luminosity L([CII]) measures the PDRs arising from 
star formation, encourages the calibration of SFR compared to L([CII]). 



Determining the SFR for dusty sources ultimately traces to the method of iKennicutt 



( 119981 ). in which the total infrared luminosity L ir is attributed to reradiation by dust of the 
primary stellar luminosity from the starbursts. The original calibration i s log SFR = log 



9.76, for L{ r in Lq. We adopt the updated calibration in Buat et a 



SFR = log L ir - 9.97. Al l of our sources have L ?r liste d in lSargsyan et al. 



(|2010f) . giving log 



( 1201 ll ) determined 



using the formulation in ISanders and Mirabel! (119961 ). described in the footnote to Table 1. 
This L- ir is an estimate of the complete 8 /iin to 1000 /mi luminosity derived from all four 
IRAS bands. 

By comparing L ir to L([CII]) for starbursts, a calibration can be determined between 
L([CII]) and SFR. It is necessary to assure that Lj r arises only from a starburst and is not 
contaminated by an AGN component. Figure 7 compares the ratio LQCII])/ L ir depending on 
source classification. The medians seen in this Figure trend as expected if L([CII]) measures 
only the starburst component but Lj r arises from both starburst and AGN components. 
As the contribution to Lj r arising from AGN luminosity increases compared to the L([CII]) 
arising from starburst luminosity, the ratio LQCII] )/Lj r systematically decreases for AGN 
(this is discussed further in section 3.3 in context of the CII deficit). For calibrating SFR to 
L([CII]) using L ir , AGN and composites in Figure 7 are not used because some fraction of 
the Li r arises from AGN. 

For starbursts only in Figure 7, the observed median ratio log L([CII])/L ir = -2.89. This 
gives the calibration log SFR = log LQCII)]) - 7.08±0.3, for SFR in M yr" 1 and LQCII]) in 
L . The uncertainty arises from the one sigma dispersion among the starburst points shown 
in Figure 7. This result is our SFR calibration. Taken with the conclusions from Figures 5 
and 6, that LQCII]) scales with the starburst component in all sources, this calibration can 
be applied to any source in which [CII] is measured. 
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Fig. 7. — Ratio of [CII] luminosity LQCII]) to total infrared luminosity L ir , compared to 
source classification from EW(PAH 6.2 /im) measured in /im. Crosses are AGN from the EW 
classification, open circles are composite AGN plus starburst, and asterisks are starbursts. 
Sources with diamonds (all AGN) are upper limits to [CII] line fluxes in Table 1. Square is 
Markarian 231. Horizontal bars are medians within each category; medians include limits 
because all limits fall below the median. Median for starbursts is used in the text to calibrate 
LQCII]) compared to star formation rate, giving log SFR = log LQCII)]) - 7.08±0.3, for SFR 
in M Q yr" 1 and LQCII]) in L . 
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The resulting measures of SFR are shown in Figure 8. The results show that the 
greatest SFRs are in sources with a starburst classification from EW(6.2 /im) even though 
these sources do not have the most luminous L ir . For example, among sources with log SFR 
> 1.8, six are starbursts, three are composite, and only two are AGN (not counting upper 
limits). This preponderance of starbursts is even greater when compared to the sample sizes; 
this high SFR includes 6 of 24 starbursts, 3 of 31 composites, and only 2 of 55 AGN. The 
dominance in SFR by starbursts arises despite the fact that the largest L ir (L ir > 10 12 L ) 
are dominated by AGN. The maximum SFR in this sample is ~ 100 M Q yr _1 , about a factor 
of ten less than in the most luminous starbursts at z ~ 2 with SFR measured using PAH 



luminosity (jWeedman and Houckll2008l ). 



Another important result from Figure 8 is that AGN are generally accompanied by 
starbursts, but AGN sources (including composites) encompass a much larger range and 
dispersion of SFR than do the pure starbursts. AGN have -1 < log SFR < 2 compared to 
0.8 < log SFR < 2.0 for starbursts. At luminosity L ir ~ 10 11 L Q , AGN have -0.2 < log 
SFR < 1.3 compared to 0.7 < log SFR < 1.7 for starbursts. Composites are intermediate. 
These results mean that some fraction of L ir arises from an accompanying starburst even 
for AGN-dominated L ir , but the large dispersion in SFK/L ir for AGN means this fraction 
varies by a factor of more than 10. 



3.3. The [CII] "Deficit" and Source Classification 



A primary result from analysis of [CII] luminosity from ISO measures was the discov- 
ery of the "[CII] Deficit " , whereby the most l uminous sources h ave weak L([CII]) compared 



to infrared lumin osity (ILuhman et al. 



Herschel results (IGracia-Carpio et al. 



20031 : iHelou et al.ll200ll ). This i s confirmed in new 



201 ll ) and ground-based results (jStacey et al.ll2010l ). 



The explanation o f this deficit remains a question, and there are many possibilities, thor- 
oughly reviewed by lLuhman et al.l (120031 ). These include HII regions with densities above the 
critical density for [CII] emission, regions with increasing ionization parameter and harder 



ioniz i ng radiation which diminishes or destroys the PDR (IMalhotra et al.ll200ll ; lAbel et al. 



2009 



Stacey et al. 201ol ; Gracia-Carpio et al. 2011 ). and self absorption or dust extinction 



of [CII]. 

One simple possibility to explain the deficit is that the most luminous sources are 
powered primarily by AGN so that most of the infrared luminosity arises from AGN which 
do not have accompanying L([CII]) from a starburst PDR. In this case, there would be 
no difference within the starburst regions or PDRs between ULIRGS and lower luminosity 
starbursts; the deficit is simply a consequence of increased AGN activity and the subsequent 
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additional L ir . This would be an important conclusion because of previous suggestions 
that individual starbursts in ULIRGS h ave higher luminosity density than other starbursts 
jLuhman et alil2003l : Istacev et allhoioh . 



Because we have already concluded that LQCII]) scales with the starburst component, 
we would expect [CII] to be relatively weaker if an AGN dominates L« r than if Li r arises 
primarily from starbursts. This is already shown in Figure 7, which demonstrates that AGN 
systematically have smaller values of L([CII])/Lj r . Whether AGN are responsible for the 
deficit is tested in Figure 9 by comparing L([CII])/ L ir with L ir using the source classifications 
derived from EW(6.2 fim). The overall dist ribution of points in Figure 9 de monstrates a 



deficit similar to the results summarized in iGracia-Carpio et al.l (j201lf ) and IStacey et al. 
( I2010l ). The line is fit to all points in our sample and would have an even steeper slope 
depending on the actual values for the upper limits. 

Figure 9 indicates that the [CII] deficit shown by our sample arises because of AGN. All 
sources with a deficit, log L([CII])/Lj r < -3.4, are either AGN or composites. All AGN do 
not show deficits, however. Some lower luminosity AGN have L([CII])/Lj r ratios similar to 
starbursts; the deficit arises only in the highest luminosity AGN. The ratios in Figure 9 and 
the SFRs shown in Figure 8 indicate that the starbursts within our sample have a maximum 
luminosity Lj r ~ 10 12 L at which log L([CII])/Lj r ~ -3.4. Any further increase in Lj r comes 
only from an AGN component without additional L([CII]), thereby decreasing L([CII])/Lj r . 

These conclusions consider only the total infrared luminosity L ir and do not address the 
important question of how the shape of the continuum spectral energy distribution depends 
on the AGN/starbur s t frac tion. A measure of "far infrared luminosity" is also defined by 
Sanders and Mirabell ( 119961 ) using L^ r as an estimate of the 40 /iin to 120 /im luminosity 
derived only from the 60 /iin and 100 /xm bands, and Lf ir is used as a luminosity measure 
in some of the other analyses cited above. Our PACS results also provide the far infrared 
continuum flux density at rest frame 158 fim, and we will discuss these results in a future 
analysis to determine, for example, if sources with larger starburst fractions determined from 
L([CII]) also have a greater proportion of far infrared luminosity. 



3.4. Comparisons to Dusty Sources at z ~ 2 

At redshifts z ~ 2, the rest frame spectral features observed with the Spitzer IRS are 
those shown in Figure 1. For starbursts and dusty, absorbed AGN, the spectra show a 
maximum near 7.8 /xm. For starbursts, this is the peak of the PAH feature. For absorbed 
AGN, this is a localized continuum maximum between absorptions on both sides of the 



=1 

CO 



O 




O 



o 



o 



o 
o 



o 



o 



o o 
oo 



o 



-o oo 

0D 



o 



o 



o 



Composite 



Starburst 



0.0 



0. 



0.2 0.3 0.4 
EW (PAH 6.2 /im 



0.5 



0.6 
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as in preceding figures. Horizontal bars are medians within each category; medians include 
limits because all limits fall below the median. 
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maximum (ISpoon et al.ll2004j ). These heavily absorbed AGN are the AGN among the DOGS 
and were generally not known from optical observations because of the heavy extinction by 
dust. Type 1 AGN or o ptically discovere d quasars are also luminous in the infrared but have 



silicate in emission (e.g. lHao et al.ll2005l ). which does not show the localized peak at 7.8 /im. 



The presence of the 7.8 jttm maximum explains why large numbers of dusty, optically 
obscured Spitzer sources are selected at z ~ 2, because this maximum becomes centered 
within the b roadband 24 ^m f ilter used for surveys with the Multiband Imaging Photometer 
for Spitzer (IRieke et al.l l20Q4f ) . Spit zer-discovered sources at these redshifts have also been 
measured in far infrared, submillimeter, and milli meter wavelengths and compared to the 



Herschel or submill imeter-discovered populations ( IMagdis et al.l I201CH : iKovacs et al.l 12010 
Lonsdale et al.l I201CH ) . (The highest redshift source confidently discovered based on IRS 
spectra has z = 3.35, a redshift limit set by the observable IRS long wavelength limit for 
identifying the 7.8 jttm maximum.) 

In spectra of faint sources with poor S/N, the spectral maximum near 7.8 iim is the 
most reliable observational measure of mid-infrared (rest frame) flux d ensity and can be 



Farrah et al. 


2008; 


Desai et al. 


2009; 


Coppin et al. 


2010) 



For these reasons, our previous 
analyses of starbursts and AGN observed with the IRS use the parameter i/L u (7.8 fim) as a 
measure of infrared luminosity for local sources and Spitzer-discovere d DOGS at z > 2 (e.g. 



Houck et al.l 120071 : IWeedman and Houckll2009bl ; ISargsvan et alJl2010h . 



Future detections or limits on [CII] at high redshifts will provide important constraints 
on the earliest epochs of formation for luminous dusty starbursts and AGN. It is useful, 
therefore, to compare f^(7.8 fim) and [CII] line fluxes for our present sample to predict [CII] 
line fluxes for the dusty population at z > 2, when scaled from the DOGS already known at 
z ~ 2. 

The comparison is sho wn in Figure 10, using the maximum flux densities f„(7.8 fiva) 
from ISargsvan et al.l (120111 ) and comparing f([CII])/z/f 1/ (7.8 /im) so that units are the same. 
If high redshift populations have the same [CII] characteristics as our low redshift Herschel 
sample, these results predict log [f([CII])/^f iy (7.8 /im)] = -2.2 ± 0.2 for starbursts and log 
[f([CII])///f v (7.8 /im)] = -2.7 ± 0.5 for AGN. The larger dispersion of this ratio for AGN 
arises because of the large range in the starburst component of AGN, discussed in section 
3.2. 
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4. Conclusions and Summary 

The Herschel PACS instrument has been used to observe [CII] 158 jum line fluxes in 
112 sources having a wide range of starburst and AGN classifications chosen because they 
have complete mid-infrared spectra with the Spitzer IRS and have complete IRAS fluxes 
for determining L ir . Of the 112 sources, 102 have reliable line detections and 10 are upper 
limits. 

It is found that the [CII] line flux correlates with the flux of the 11.3 ^m PAH feature, log 
[f([CII] 158 /im)/f(11.3 /im PAH)] = -0.22 ± 0.25. This f([CII])/f(PAH) ratio is independent 
of AGN/starburst classification as determined from equivalent width of the 6.2 /xm PAH 
feature. We conclude that [CII] line flux measures the starburst component of any source as 
reliably as the PAH feature. 

This conclusion leads to a calibration of star formation rate determined from the lumi- 
nosity of [CII] for the starburst component in any source . The calibration is derived using 
L ir only for starbursts to avoid AGN contamination of L ir and has the result log SFR = log 
LQCII)]) - 7.08 ± 0.3, for SFR in M yr" 1 and LQCII]) in L . This result applies to the 
starburst component of any source in which [CII] is observed. The maximum SFRs in the 
sample are 100 M Q yr" 1 , and SFRs are dominated by sources classified as starbursts, but 
most AGN also have some measurable starburst component. 

The [CII] "deficit", or a smaller ratio of L([CII])/Lj r with increasing Lj r , is shown to 
arise because L ir of the most luminous sources arises primarily from an AGN so that LQCII]) 
from the starburst component is small in comparison. 
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Table I. [CII] Line Fluxes and Luminosities 



No. 


l\Tn m r> 


rrirwfi l n attic 






EW a 
6 . 2 fim 




S/N c 


fCTTl /1 1 3 d 


fCTTl n 8 e 


^CII 


^CII 1 ^ir 


1J.CI 






J2000 








W m~ 2 




iU e> 


log 


log 


log 




1 


Mrk0334 


000309.62+215736.6 





.0219 


0.248 


6.70c-16 


169 


-0.26 


-2.27 


8.27 


-2.81 


1342235847 


2 


MCG-02-01-051/2 


001850.90-102236.7 


0. 


.0271 


0.558 


1.18c-15 


4.9 


-0.17 


-2.09 


8.70 


-2.75 


1342235846 


3 


IRAS00199-7426 


002207.01-740941.7 





.0964 


0.351 


1.98c-16 


51 


-0.29 


-2.30 


9.05 


-3.27 


1342237574 


4 


E12-G21 


004046.10-791424.0 





0330 


0.278 


5.76c-16 


5.7 


-0.12 1 


-2.04 


8.48 


-2.54 


1342232583 


5 


IRASF00456-2904SW 


004806.75-284818.6 





.1103 


0.428 


1.26c-16 


44 


-0.25 


-2.19 


8.98 


-3.23 


1342238141 


6 


MCG-03-04-014 


011008.93-165109.9 





.0350 


0.455 


1.38c-15 


5.2 


-0.25 


-2.14 


8.95 


-2.67 


1342238385 


7 


NGC0454 


011424.90-552352.0 





.0122 


0.008 k 


1.55e-16 


5.2 


-0.31 


-2.51 


7.11 


-3.14 


1342232616 


8 


ESO244-G012 


011808.31-442743.4 





.0229 


0.536 


1.14c-15 


188 


-0.32 


-2.26 


8.53 


-2.86 


1342234998 


9 


ESO353-G020 


013451.26-360814.4 


0. 


0159 


0.404 


1.56e-15 


249 


-0.18 


-2.19 


8.35 


-2.68 


1342238600 


10 


IRASF01364-1042 


013852.91-102711.0 





.0482 


0.306 


1.32c- 16 


6.0 


-0.05 


-2.27 


8.26 


-3.54 


1342238601 


11 


UGC01385 


015453.82+365504.3 





.0188 


0.538 


5.70c-16 


146 


-0.44 


-2.29 


8.06 


-2.97 


1342237475 


12 


NGC0788 


020106.45-064855.9 





.0136 


<0.003 k 


8.43c-17 


4.8 


0.12 


-2.56 


6.95 


-3.18 


1342238364 


13 


IRAS02054+0835 j 


020806.90+085004.3 





.3450 


0.024 k 


<2.00c-17 


0.1 


<-0.30 


<-2.88 


<9.26 


<-3.80 


1342239376 


14 


Mrk0590i 


021433.56-004600.1 





.0264 


<0.007 


2.06c-16 


3.5 


0.84 


-1.86 


7.92 


-2.75 


1342238636 


15 


UGC01845 


022407.97+475811.9 





.0156 


0.443 


1.42e-15 i 


3.4 


-0.37 


-2.32 


8.29 


-2.80 


1342239504 


16 


IC1816 


023151.00-364019.4 





.0169 


0.031 k 


3.42c- 16 


52 


0.12 


-1.89 


7.75 


-2.72 


1342239370 


17 


NGC0973 


023420.11+323020.2 





.0162 


<0.002 k 


2.46C-16 1 


3.9 


0.47 


-2.05 


7.56 


-2.89 


1342239500 


18 


IRASF02437+2122 


024639.13+213510.4 





.0233 


0.155 


1.50e-16 


44 


-0.42 


-2.50 


7.67 


-3.49 


1342239499 


19 


UGC02369 


025401.84+145815.7 





.0312 


0.434 


8.21C-16 1 


174 h 


-0.10 


-2.08 


8.66 


-2.96 


1342239497 


20 


Mrkl066 


025958.59+364914.3 





.0120 


0.249 


9.39e-16 


4.9 


-0.46 


-2.30 


7.89 


-3.04 


1342238915 


21 


IRASF03217+4022 


032505.37+403332.2 





.0234 


0.406 


5.96C-16 1 


3.7 


-0.29 


-2.21 


8.27 


-3.04 


1342238940 


22 


Mrk0609 


032525.34-060838.7 





.0345 


0.318 


5.44e-16 


4.0 


-0.31 


-2.12 


8.57 


-2.81 


1342239752 


23 


IRASF03359+1523 


033847.07+153254.1 





.0354 


0.358 


4.60c- 16 


5.0 


-0.04 


-1.93 


8.52 


-3.04 


1342238916 


24 


IRASF03450+0055 


034740.18+010514.0 





.0310 


<0.001 


2.69c-17 


6.8 


>-0.42 


-3.33 


7.17 


-3.87 


1342238943 


25 


IRAS03538-6432 


035425.23-642344.5 





.3007 


0.079 k 


<2.00c-17 


2.4 


<-0.55 


<-2.73 


<9.13 


<-3.64 


1342223119 


26 


IRAS04103-2838 


041219.53-283024.4 





.1175 


0.182 


8.53e-17 


23.5 


-0.41 


-2.46 


8.87 


-3.33 


1342239509 


27 


IRAS04114-5117 


041244.92-510934.2 





.1250 


0.338 


6.32c-17 


4.4 h 


-0.25 


-2.29 


8.79 


-3.46 


1342226902 


28 


ESO420-G013 


041349.70-320025.3 





.0119 


0.208 


1.38c-15 


4.0 


-0.32 


-2.29 


8.04 


-3.00 


1342238379 


29 


3C120 


043311.10+052115.6 





.0330 


<0.001 


3.62e-16 


6.3 


0.73 


-2.27 


8.36 


-2.82 


1342225795 


30 


ESO203-IG001J 


044649.55-483330.6 





.0529 


<0.029 k 


5.14c-17 


5.1 


0.06 


-2.66 


7.93 


-3.91 


1342238378 


31 


MCG-05-12-006 


045204.96-325926.0 





.0188 


0.400 


3.74c-16 


107 


-0.59 


-2.42 


7.88 


-3.28 


1342239733 


32 


Arkl20 


051611.42-000859.4 





.0327 


<0.001 


1.34e-16 


6.2 


0.00 


-2.78 


7.92 


-3.09 


1342226750 


33 


VIIZw31 


051646.39+794012.9 





.0539 


0.491 


8.55e-16 


4.0 


-0.33 


-2.23 


9.16 


-2.81 


1342219853 



Table 1 — Continued 



No. Name coordinates z EW a f[CII] b S/N c [CII]/11.3 d [CII]/7.8 C L C n C L C u/L ir ^ Herschel id 







J2000 






6.2 (im 
fim 


w m - 2 




log 


log 


solar 
log 


log 




34 


IRASF05187-1017 


052106.53-101446.2 





.0283 


0.414 


2.20c-16 


6.0 


-0.25 


-2.25 


8.01 


-3.26 


1342227348 


35 


2MASXJ05580206-3820043i 


055802.00-382004.0 


0. 


0339 


<0.001 k 


1.29c-17 


3.8 


>-0.44 


-4.05 


6.93 


-4.23 


1342239744 


36 


IRASF06076-2139 


060945.74-214024.5 


0. 


.0374 


0.238 


l.llc-16 


4.2 


-0.51 


-2.61 


7.96 


-3.67 


1342230910 


37 


IRAS06301-7934-> 


062642.20-793630.4 





.1564 


0.148 


3.10c-17 


3.9 


-0.10 


-2.39 


8.69 


-3.68 


1342231278 


38 


IRAS06361-6217j 


063635.71-622031.8 


0, 


.1596 


0.052 k 


3.38e-17 


3.6 


-0.09 


-2.80 


8.75 


-3.66 


1342238377 


39 


NGC2273 


065008.72+605045.0 





0061 


0.126 


7. 18c- 16 


161 


-0.44 


-2.42 


7.18 


-2.98 


1342230996 


40 


UGC03608 


065734.41+462410.6 





.0214 


0.424 


1.46c-15 


212 


0.27 


-1.64 


8.46 


-2.72 


1342230955 


41 


IRASF06592-6313 


065940.26-631752.4 





0230 


0.402 


2.51c-16 


5.5 


-0.62 


-2.47 


7.88 


-3.31 


1342231286 


42 


Mrk0009 


073657.00+584613.0 


0. 


.0399 


<0.09' 


1.12c- 16 


4.2 


-0.50 


-2.73 


8.02 


-3.11 


1342243533 


43 


IRAS07598+6508J 


080430.46+645952.9 





.1488 


0.005 


2.87c-17 


3.9 


-0.02 


-3.40 


8.61 


-3.90 


1342243534 


44 


Mrk0622 


080741.04+390015.2 





.0232 


0.245 


8.36c-17 


4.3 


-0.63 


-2.33 


7.41 


-3.31 


1342229688 


45 


ESO60-IG016 


085232.07-690154.8 





.0455 


0.070 


4.99e-16 


113 


0.18 


-2.25 


8.77 


-2.95 


1342228521 


46 


Mrk0018 


090158.39+600906.0 





.0111 


0.275 


6.47c- 16 


5.1 


-0.10 


-1.93 


7.65 


-2.52 


1342231958 


47 


MCG-01-24-012 


092046.25-080322.1 





0196 


<0.005 k 


1.39c- 16 


5.2 


0.30 


-2.26 


7.49 


-2.87 


1342231718 


48 


Mrk0705 


092603.29+124403.6 


0. 


0292 


0.018 


1.34e-16 


5.1 


0.00 


-2.37 


7.82 


-2.88 


1342231715 


49 


IRASF10038-3338 


100604.65-335306.1 


0. 


.0342 


<0.002 k 


4.03e-16 


86 


0.05 


-2.89 


8.44 


-3.25 


1342235705 


50 


NGC3393 


104823.39-250942.8 





.0125 


<0.003 


3.44e-16 


70 


0.73 


-1.88 


7.48 


-2.91 


1342232587 


51 


IRAS11119+3257 


111438.88+324133.1 





.1890 


0.014 k 


<2.00e-17 


3.0 


<-0.62 


<-3.30 


<8.68 


<-3.94 


1342232307 


52 


ESO319-G022 


112754.18-413651.7 





.0164 


0.302 


4.00e- 16 


72 


-0.05 


-2.00 


7.78 


-3.22 


1342235704 


53 


IRAS12018+194lj 


120424.53+192509.8 





.1686 


0.108 


2.41c-17 


7.9 


-0.11 


-2.79 


8.65 


-3.86 


1342223403 


54 


UGC07064 


120443.34+311038.2 


0. 


.0250 


0.081 


3.94e-16 


4.1 


0.10 


-1.91 


8.15 


-2.89 


1342223404 


55 


NGC4507 


123536.55-395433.3 


0. 


.0118 


<0.002 


5.62e-16 


5.7 


0.09 


-2.46 


7.65 


-2.98 


1342225720 


56 


PG1244+026J 


124635.24+022208.7 





.0482 


0.017 


1.65c-17 


3.5 


-0.27 


-2.80 


7.35 


-3.65 


1342236279 


57 


IRAS12514+1027 j 


125400.82+101112.4 





.3182 


<0.009 k 


<2.00c-17 


1.6 


<-0.29 


<-3.00 


<9.18 


<-3.60 


1342237582 


58 


ESO507-G070 


130252.42-235517.8 





.0217 


0.417 


9.51e-16 


165 


0.11 


-2.07 


8.41 


-3.03 


1342225742 


59 


NGC4941 


130413.10-053306.0 





.0037 


<0.041 k 


1.71c-16 


4.9 


0.32 


-2.09 


6.12 


-3.26 


1342225749 


60 


ESO323-G077 


130626.13-402452.0 





.0156 


0.029 


8.27c-16 


5.4 


-0.17 


-2.41 


8.02 


-2.93 


1342225811 


61 


MCG-03-34-064 


132224.45-164342.4 


0. 


.0165 


<0.001 k 


1.90e-16 


6.2 


-0.28 


-3.04 


7.47 


-3.65 


1342225746 


62 


IRASF13279+3401 j 


133015.23+334629.4 


0. 


.0238 


0.038 k 


2.58e-17 


3.4 


-0.13 


-2.55 


6.83 


-3.63 


1342232550 


63 


M-6-30-15 


133554.54-341750.4 





.0077 


<0.001 


7.04e-17 


4.4 h 


-0.49 


-3.11 


6.38 


-3.53 


1342225810 


64 


IRAS13342+3932 


133624.07+391730.1 


0, 


.1793 


0.073 


6.54c-17 


5.2 


-0.20 


-2.49 


9.15 


-3.30 


1342234948 


65 


IRAS13352+6402J 


133651.15+634704.7 





2366 


0.076 k 


<2.00c-17 


2.4 


<-0.40 


<-2.80 


<8.89 


<-3.67 


1342231427 


66 


IRASF13349+2438i 


133718.73+242303.3 





.1076 


<0.001 k 


2.54e-17 


4.9 


-0.24 


-3.79 


8.26 


-4.04 


1342236984 



Table 1 — Continued 



No. Name coordinates z EW a f[CII] b S/N c [CII]/11.3 d [CII]/7.8 e Leu 1 L C u/L ir ^ Herschel id 







J2000 






6.2 (im 
fim 


W m~ 2 




log 


log 


solar 
log 


log 




67 


NGC5347 


135317.83+332927.0 





.0078 


0.025 


1.15c- 16 


29.0 


-0.40 


-2.66 


6.60 


-3.38 


1342223133 


68 


IRAS14026+4341J 


140438.72+432707.3 


0. 


3233 


0.011 k 


<2.00c-17 


0.1 


<-0.33 


<-2.92 


<9.20 


<-3.71 


1342236271 


69 


OQ+208 


140700.39+282714.0 


0. 


.0766 


0.018 


6.28c-17 


7.1 


-0.43 


-2.74 


8.35 


-3.32 


1342237584 


70 


NGC5548 


141759.53+250812.4 





.0172 


0.024 


1.70e-16 


4.5 


-0.12 


-2.48 


7.45 


-3.17 


1342236985 


71 


Mrkl490 


141943.27+491411.9 


0. 


.0257 


0.448 


3.18e-16 


103 


-0.60 


-2.48 


8.08 


-3.22 


1342232331 


72 


PG1426+015 j 


142906.59+011706.5 


0. 


.0865 


<0.002 


2.30e-17 


3.2 


-0.14 


-3.06 


8.02 


-3.48 


1342238159 


73 


PG1440+356 


144207.46+352622.9 





.0780 


0.017 


5.77c-17 


13.2 


-0.37 


-2.69 


8.34 


-3.19 


1342223736 


74 


NGC5728 


144223.93-171511.0 





.0094 


0.077 


1.25e-15 


6.4 


0.07 


-1.85 


7.79 


-2.80 


1342238133 


75 


NGC5793 


145924.76-164136.0 


0. 


.0116 


0.393 


8.98c- 16 1 


3.6 


-0.21 


-2.24 


7.84 


-2.79 


1342238134 


76 


IRAS15001+1433 


150231.94+142135.3 





1627 


0.140 


6.40c-17 


5.6 


-0.32 


-2.50 


9.05 


-3.42 


1342236887 


77 


IRAS15225+2350 j 


152443.94+234010.2 





.1390 


0.051 k 


3.23e-17 


3.2 


-0.17 


-2.65 


8.60 


-3.57 


1342238152 


78 


Mrk0876 


161357.18+654309.6 


0. 


.1290 


0.010 


2.99c-17 


3.6 


-0.51 


-2.98 


8.50 


-3.47 


1342222163 


79 


IRASF16164-0746 


161911.75-075403.0 





0235 


0.439 


7.51C-16 1 


4.6 


-0.11 


-2.23 


8.50 


-3.06 


1342227799 


80 


Mrk0883 


162952.85+242638.3 





.0378 


0.263 


1.84c- 16 


5.5 


-0.10 


-1.83 


8.18 


-2.77 


1342238911 


81 


CGCG052-037 


163056.53+040458.7 


0. 


.0245 


0.510 


1.13c-15 


233 


-0.20 


-2.12 


8.59 


-2.78 


1342238909 


82 


IRAS16334+4630 


163452.37+462453.0 


0. 


1908 


0.302 


3.62c-17 


4.2 


-0.47 


-2.44 


8.95 


-3.50 


1342232265 


83 


ESO069-IG006 


163811.85-682608.2 





.0464 


0.497 


1.080-15' 


3.1 


-0.27 


-2.27 


9.14 


-2.79 


1342231301 


84 


IRASF16399-0937N 


164240.11-094313.7 


0. 


.0270 


0.175 


1.49e-15 


40 


0.39 


-1.62 


8.80 


-2.70 


1342240777 


85 


2MASSJ165939.77+183436.9 j 


165939.77+183436.9 


0. 


.1707 


0.016 


<2.00c-17 


2.4 


<-0.30 


<-2.95 


<8.58 


<-3.56 


1342238910 


86 


PG1700+518 


170124.91+514920.4 





.2920 


<0.009 


<2.00c-17 


0.1 


<-0.56 


<-3.11 


<9.10 


<-3.61 


1342225993 


87 


IRAS17044+6720J 


170428.41+671628.5 





.1350 


0.040 k 


2.06c-17 


3.3 


-0.39 


-3.07 


8.38 


-3.79 


1342223716 


88 


IRAS17068+4027 


170832.12+402328.2 


0. 


1790 


<0.079 k 


3.37c-17 


3.4 


-0.32 


-2.68 


8.86 


-3.54 


1342232574 


89 


IRASF17132+5313 


171420.45+531031.6 





.0509 


0.420 


6.090-16' 


3.9 h 


-0.12 


-2.00 


8.97 


-2.91 


1342223740 


90 


ESO138-G027 


172643.35-595555.2 


0. 


.0208 


0.351 


1.44e-15 


3.9 


0.19 


-1.76 


8.55 


-2.78 


1342240168 


91 


CGCG141-034 


175656.65+240102.0 


0. 


0198 


0.388 


6.65c-16 


7.8 


-0.17 


-2.18 


8.17 


-2.93 


1342231759 


92 


H1821+643 


182157.31+642036.3 





.2970 


<0.002 


<2.00c-17 


2.5 


<-0.57 


<-3.37 


<9.11 


<-4.03 


1342222100 


93 


IC4734 


183825.75-572925.4 





0156 


0.380 


1.36e-15 


4.8 


-0.12 


-2.10 


8.28 


-2.99 


1342241710 


94 


IRAS18443+7433 j 


184254.80+743621.0 


0. 


1347 


<0.051 k 


2.31e-17 


3.7 


-0.21 


-2.92 


8.43 


-3.89 


1342232253 


95 


ESO140-G043 


184453.98-622153.4 


0. 


0142 


0.022 


1.76c-16 


45 


-0.59 


-2.80 


7.30 


-3.29 


1342240169 


96 


lH1836-786 j 


184703.20-783151.0 


0. 


0741 


<0.012 k 


2.49c-17 


3.4 


-0.20 


-2.75 


7.92 


-3.49 


1342231316 


97 


ESO593-IG008 


191431.15-211906.3 


0. 


0487 


0.405 


9.51C-16 1 


3.1 


0.02 


-2.02 


9.12 


-2.76 


1342231749 


98 


ESO-141-G055 


192114.15-584013.1 


0. 


0371 


<0.003 


1.32c-16 


4.3 


0.24 


-2.54 


8.02 


-3.04 


1342231317 


99 


ESO339-G011 


195737.60-375608.4 


0. 


.0192 


0.188 


1.54e-15 


5.8 


0.22 


-1.87 


8.51 


-2.62 


1342232296 



Table 1 — Continued 



No. Name coordinates z EW a f[CII] b S/N c [CII]/11.3 d [CII]/7.8 C L C u C L C n/L ir s Herschel id 

6.2 fim solar 







J2000 








W m" 2 




log 


log 


log 


log 




100 


IRAS20037-1547 


200631.70-153908.0 


0. 


1919 


0.043 


4.63e-17 


3.3 


-0.32 


-2.74 


9.06 


-3.56 


1342232297 


101 


NGC6860 


200846.90-610601.0 


0. 


.0149 


0.084 


3.69e-16 


6.3 


0.01 1 


-2.22 


7.67 


-2.67 


1342232295 


102 


ESO286-G035 


210411.11-433536.1 


0. 


0174 


0.504 


1.63e-15 


4.7 


-0.19 


-2.13 


8.45 


-2.63 


1342232563 


103 


NGC7213 


220916.25-471000.0 


0. 


0058 


<0.005 


2.74e-16 


5.3 


0.03 


-2.22 


6.72 


-3.25 


1342232569 


104 


ESO602-G025 


223125.48-190204.0 


0. 


0250 


0.300 


1.24e-15 


190 


-0.19 


-2.18 


8.65 


-2.69 


1342233481 


105 


UGC12138 


224017.05+080314.1 


0. 


.0250 


0.081 


1.43e-16 


4.1 


-0.43 


-2.37 


7.71 


-3.06 


1342235675 


106 


UGC12150 


224112.21+341456.8 


0. 


.0214 


0.383 


8.96e-16 


5.6 


-0.29 


-2.20 


8.37 


-2.95 


1342223714 


107 


ESO239-IG002 


224939.84-485058.3 


0. 


0430 


0.309 


1.91e-16 


6.2 


-0.50 


-2.44 


8.31 


-3.49 


1342232568 


108 


Zw453.062 


230456.55+193307.1 


0. 


.0251 


0.379 


1.05e-15 


5.3 


0.18 


-1.82 


8.58 


-2.78 


1342235676 


109 


IRAS23060+0505 


230833.97+052129.8 


0. 


1730 


0.009 k 


<2.00e-17 


3.0 


<-0.48 


<-3.49 


<8.60 


<-3.89 


1342235673 


110 


NGC7603 


231856.62+001438.2 


0. 


.0295 


0.018 


2.63e-16 i 


4.2 


-0.28 


-2.62 


8.12 


-2.72 


1342222576 


111 


MCG-83-1 


234200.91-033654.4 


0. 


0232 


0.504 


1.47e-15 


5.8 


0.09 


-1.85 


8.66 


-2.52 


1342235844 


112 


CGCG381-051 


234841.70+021423.0 


0. 


.0307 


0.542 


2.34e-16 


5.1 


-0.24 1 


-2.13 


8.10 


-3.07 


1342236876 



a Rest frame equivalent width of 6.2 /im PAH emission feature from lSaresvan et al. (2011) used to classify source as AGN, composite, or starburst. I 

b Total flux of [CII] 158 fim emission line in units of W m~ 2 using gaussian fit to line for simple profiles and integrated flux of line for complex profiles 
(noted by footnote i in this column). Line flux listed is the total flux observed within the 9 spaxels centered on the brightest spaxel, increased by a correction 
factor of 1.16 to 1.21 to include flux that would fall outside these spaxels for an unresolved source. The correction factor adopted for the range of observed 
[CII] wavelengths from 160 fim to 210 fim is 1.16(A/158 /im) 0,17 . Uncertainties of individual fits given by S/N in next column; systematic uncertainty for 
all fluxes depends on PACS flux calibration, estimated as ± 12% in the PACS Spectroscopy performance and calibration document PICC-KL-TN-041 . 

c Signal to noise ratio of total line flux in brightest spaxel, using one sigma uncertainty of profile fit. 

d Ratio of flux in [CII] 158 fim to flux of PAH 11.3 fim emission feature, from lSaresvan et al. I ll201ll ). 

e Ratio f([CII])/i/f„(7.8 /im) using 1^(7.8 fim) from lSaresvan et al.l iboill) . 

[C II] 158 fim emission line luminosity in Lq using luminosity distances determined for Ho = 71 kms" 1 Mpc -1 , Sljv/ =0. 27 and £2^=0.73, from IWrightl 
d2006T l : |http:/ /www. astro. ucla.edu/~wright/CosmoCalc. html| [log Lqii (solar) = log Lqjj (W) - 26.59]. 

g Ratio of [CII] luminosity to Li r using Li r given in ISargsvan et al.l ll201ll ) from f^ r determined as in I Sanders and Mirabell dl996T) . f; r = 1.8 x 
10" 11 [13.48f, y (12) + 5.16f„(25) + 2.58f„(60) + f„(100)] in erg cm" 2 s" 1 using IRAS flux densities at 12 fim, 25 fim, 60 /im and 100 /mi. 

^Brightest spaxel displaced one spaxel from central 3,3 spaxel. 

'[CII] line profile is asymmetric or has component structure so total line flux is integrated flux including all components rather than flux within a single 



gaussian fit. 

j Observation made with 2 repetition cycles of line spectroscopy point source chop nod mode; all observations without note made with single cycle. 
k IRS spectrum shows 9.7 /an silicate feature in absorption. 

'For these 4 sources, spectra are not in CASSIS because IRS spectra obtained in mapping mode; PAH measures from lWu et al. 
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